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Abstract 

The exclusive semileptonic B — ► n(p)ev decays are studied in the framework of the rela- 
tivistic quark model based on the quasipotential approach in quantum field theory. The 
large recoil momentum of the final 7r(p) meson allows for the expansion of the decay form 
factors at q 2 = in the inverse powers of 6-quark mass. This considerably simplifies 
the analysis of these decays. The l/m& expansion is carried out up to the first order. 
It is found that the g 2 -dependence of the axial form factor A\ is different from the q 2 - 
behaviour of other form factors. We find T(B — > ixev) = (3.1 ±0.6) x \V ub \ 2 x 10 12 s _1 and 
T(B — > peu) = (5.7 ± 1.2) x \V u b\ 2 x 10 12 s _1 . The relation between semileptonic and rare 
^ ; radiative £>-decays is discussed. 
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1 Introduction 
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The investigation of semileptonic decays of B mesons into light mesons is important for the 
determination of the Cabibbo-Kobayashi-Maskawa matrix element V u b, which is the most 
poorly studied. At present the value of V u b is mainly determined from the endpoint of the 
lepton spectrum in smileptonic .B-decays [1] . Unfortunately, the theoretical interpretation 
of the endpoint region of the lepton spectrum in inclusive B — > X u iv decays is very 
complicated and suffers from large uncertainties [2]. The other way to determine V u \> 
is to consider exclusive semileptonic decays B — > ii(p)ev. These are the heavy-to-light 
transitions with a wide kinematic range. In contrast to the heavy-to-heavy transitions, 
here we can not expand matrix elements in the inverse powers of the final quark mass. 
It is also necessary to mention that the final meson has a large recoil momentum almost 
in the whole kinematical range. Thus the motion of final 7r(p) meson should be treated 
relativistically. If we consider the point of maximum recoil of the final meson, we find 
that 7r(p) bears the large relativistic recoil momentum |A maa .| of order of m&/ '2 and the 
energy of the same order. Thus at this kinematical point it is possible to expand the 
matrix element of the weak current both in inverse powers of 6-quark mass of the initial 
B meson and in inverse powers of the recoil momentum |A ma:r | of the final 7r(p) meson. 
As a result the expansion in powers l/m& arises for the B — > ir(p) semileptonic form factor 
at q 2 = 0, where q 2 is a momentum carried by the lepton pair. The aim of this paper 
is to realize such expansion in the framework of relativistic quark model. We show that 
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this expansion considerably simplifies the analysis of exclusive B — > it(p)ev semileptonic 
decays. 

Our relativistic quark model is based on the quasipotential approach in quantum field 
theory with the specific choice of the qq potential. It provides a consistent scheme for 
calculation of all relativistic corrections at a given order of v 2 /c 2 and allows for the heavy 
quark 1 / rriQ expansion. This model has been applied for the calculations of meson mass 
spectra [3], radiative decay widths [4], pseudoscalar decay constants [5], heavy-to-heavy 
semileptonic [6] and nonleptonic [7] decay rates. The heavy quark \jvriQ expansion in 
our model for the heavy-to-heavy semileptonic transitions has been developed in [8] up 
to l/rriQ order. The results are in agreement with the model independent predictions 
of the heavy quark effective theory (HQET) [9]. The \jm b expansion of rare radiative 
decay form factors of B mesons has been carried out in [10] along the same lines as in the 
present paper. 



2 Relativistic quark model 

In the quasipotential approach meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [11] of the Schrodinger 
type [12]: 

where the relativistic reduced mass is 

!>* ~ up > ( 2 ) 

u2< A/r\ _ i M " ~ ( m - + m b f][M 2 - (m a - m b f] 

b [M) ~ 4M^ ' [6) 

m a ,b are the quark masses; M is the meson mass; p is the relative momentum of quarks. 
While constructing the kernel of this equation V(p, q; M) — the quasipotential of quark- 
antiquark interaction — we have assumed that effective interaction is the sum of the 
one-gluon exchange term with the mixture of long-range vector and scalar linear confin- 
ing potentials. We have also assumed that at large distances quarks acquire universal 
nonperturbative anomalous chromomagnetic moments and thus the vector long-range po- 
tential contains the Pauli interaction. The quasipotential is defined by [3]: 

V(p, q, M) = u a {p)u b {-p){^asD^)^ al l^ (4) 

where as is the QCD coupling constant, is the gluon propagator; 7^ and u{p) are the 
Dirac matrices and spinors; k = p — q; the effective long-range vector vertex is 

r„(k) = 7 M + ( 5 ) 
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k is the anomalous chromomagnetic quark moment. Vector and scalar confining potentials 
in the nonrelativistic limit reduce to 

V c v onf (r) = (l-e)(Ar + B), V? onf {r) = e{Ar + 5), (6) 

reproducing V^°"/ d (r) = V^ on ^ + V^ on ^ = Ar + B, where e is the mixing coefficient. The 
explicit expression for the quasipotential with the account of the relativistic corrections 
of order v 2 /c 2 can be found in ref. [3]. All the parameters of our model: quark masses, 
parameters of linear confining potential A and B, mixing coefficient e and anomalous 
chromomagnetic quark moment k were fixed from the analysis of meson masses [3] and 
radiative decays [4]. Quark masses: m& = 4.88 GeV; m c = 1.55 GeV; m s = 0.50 GeV; 
m u,d — 0.33 GeV and parameters of linear potential: A = 0.18 GeV 2 ; B = —0.30 GeV 
have standard values for quark models. The value of the mixing coefficient of vector 
and scalar confining potentials e = —0.9 has been primarily chosen from the considera- 
tion of meson radiative decays, which are very sensitive to the Lorentz-structure of the 
confining potential: the resulting leading relativistic corrections coming from vector and 
scalar potentials have opposite signs for the radiative Ml-decays [4] . Universal anomalous 
chromomagnetic moment of quark k — — 1 has been fixed from the analysis of the fine 
splitting of heavy quarkonia 3 Pj- states [3]. 

Recently we have considered the expansion of the matrix elements of weak heavy 
quark currents between pseudoscalar and vector meson states up to the second order in 
inverse powers of the heavy quark masses [8] . It has been found that the general structure 
of leading, subleading and second order 1/toq corrections in our relativistic model is in 
accord with the predictions of HQET. The heavy quark symmetry and QCD impose rigid 
constraints on the parameters of the long-range potential of our model. The analysis of 
the first order corrections [8] allowed to fix the value of effective long-range anomalous 
chromomagnetic moment of quarks k — — 1, which coincides with the result, obtained from 
the mass spectra [3] . The mixing parameter of vector and scalar confining potentials has 
been found from the comparison of the second order corrections to be e — — 1. This value 
is very close to the previous one e = —0.9 determined from radiative decays of mesons 
[4]. Therefore, we have got QCD and heavy quark symmetry motivation for the choice 
of the main parameters of our model. The found values of e and k imply that confining 
quark- ant iquark potential has predominantly Lorentz-vector structure, while the scalar 
potential is anticonfining and helps to reproduce the initial nonrelativistic potential. 

3 B — > Ti{p)ep decay form factors 

The form factors of the semileptonic decays B — > -neu and B — > peu are defined in the 
standard way as: 

{^{P^)\qiM B (PB)) = f+(Q 2 )(PB+Pn)^ + f-(q 2 )(PB-Pir)n, (7) 

(p( Pp , e)\q^(l - i 5 )b\B(p B )) = -(M B + MjAtfty 
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' M B W^^^^W (8) 
where g = PB~Pn(p), e is a polarization vector of p meson. In the limit of vanishing lepton 
mass, the form factors /_ and A 3 do not contribute to the decay rates and thus will not 
be considered. 

The matrix element of the local current J between bound states in the quasipotential 
method has the form [13]: 

(n(p)\J,(0)\B) = J ^^^(p)(p)r,( P ,q)^(q), (9) 

where T At (p, q) is the two-particle vertex function and Vf^s are the meson wave functions 
projected onto the positive energy states of quarks. 

In the case of semileptonic decays J M = 97 M (1 — 7 5 )6 and in order to calculate its 
matrix element between meson states it is necessary to consider the contributions to T 
from Figs. 1 and 2. Thus the vertex functions look like 

r«(p, q) = u,(pi) 7 „(l - 7 5 H(<7i)(27r) 3 5(p 2 - q 2 ), (10) 

and 



lf(p, q) = ^( Pl K(P2) {71,(1 - 7i 5 ) ^X^ 7 1 V( P2 - q 2 ) 



£ b {ki) + e b {pi) 

~ (2T!\ s ^^ 1 ~ ^WuhM, (ii) 

e q {Ki) + £q\qi) 

where ki = px - A; k' x = q x + A; A = p B -p 7r(p) ; e(p) = (m 2 + p 2 ) 1/2 ; 

A(")(p) 



e(p) - (m7° + 7°(7p)) 



2e(p) 

Note that the contribution I^ 2 ) is the consequence of the projection onto the positive- 
energy states. The form of the relativistic corrections resulting from the vertex function 
I^ 2 ) is explicitly dependent on the Lorentz-structure of gg-interaction. 

It is convenient to consider the decay B — > -n(p)ev in the B meson rest frame. Then the 
wave function of the final 7r(p) meson moving with the recoil momentum A is connected 
with the wave function at rest by the transformation [13] 

^(p)a(p) = <Wa)<Wa)^(p)o(p), (12) 

where D 1 ^ 2 (R) is the well-known rotation matrix and R w is the Wigner rotation. 

The meson wave functions in the rest frame have been calculated by numerical solution 
of the quasipotential equation (1) [14]. However, it is more convenient to use analytical 
expressions for meson wave functions. The examination of numerical results for the ground 
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state wave functions of mesons containing at least one light quark has shown that they 
can be well approximated by the Gaussian functions 




with the deviation less than 5%. 
The parameters are 



p B = 0.41 GeV; (3 <p) = 0.31 GeV. 

Substituting the vertex functions (10), (11), with the account of the Lorentz trans- 
formation of final meson wave function (12), in the matrix element (9) we can get the 
expressions for the semileptonic decay form factors defined in eqs. (7) and (8). They are 
very cumbersome and lengthy especially those which arise from the vertex function 
in Fig. 2, because they explicitly depend on the form of the gg-interaction potential and 
it is necessary to integrate both with respect to d 3 p and d 3 q (see eq. (9)). However, as it 
was already mentioned in the introduction, at the point of maximum recoil of final 7r(p) 
meson the large value of recoil momentum |A ma:E | ~ ?7i&/2 allows for the expansion in 
powers of l/m b for the decay form factors /+(0), ^i(O), A 2 (0) and V(0). This expansion 
leads to considerable simplification of the formulae. The large value of recoil momentum 
I A. max \ permits to neglect p 2 in comparison with in the quark energy e q (p + A) of 

final meson in the expressions for the form factors originating from []. Thus we can 
perform one of the integrations in the current matrix element (9) using the quasipoten- 
tial equation as in the case of the heavy final meson [4,6]. As a result, we get compact 
formulae suitable for numerical treatment. Neglecting terms of the second order in l/ra& 
expansion we get for the form factor of B — > 7rez/ decay: 





E n + M w max J\ 2e q (p + A 



e q (p + A max ) +m q 



, 3 ( P A ma J | pi+p 2 y 



s q (p + A max )+m q 2 A? max (e q (p) + m q ) (e q (p + A max ) + m q ) 
pI+pI 2 (pA mo J 1 / M B -E n \ 



1 + 



e q (p) + m q ) 2 £ b (A max ) + m b A? max e q (p) + m q \ e b (A max ) + m h 



x (m b + M n - e b (p) - e q {p) - 2e q [p + - ^ A^)) jtt B (p), (14) 

where we have set [3] the mixing parameter (6) of vector and scalar confining potentials 
e = —1 and long-range anomalous chromomagnetic quark moment (5) k = —1; 

|A_| = M °- M -; E, = M * + M *; (15) 

I max , n , V ) 



1 Such approximation corresponds to the omitting terms of the third order in 1/ m,b expansion. 
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and z-axis is chosen in the direction of A. 

Similar expressions can be written for the form factors of B — > pev decay at q 2 = 0. 
They have the same structure as (14) and will be given elsewhere [15]. 

Let us proceed further and for the sake of consistensy carry out the complete expansion 
of form factor (14) in inverse powers of 6-quark mass. 

The mass of B meson has the following expansion in l/m b [9] 

M B = m b + A + o(—), (16) 

In our model parameter A is equal to the mean value of the light quark energy inside the 
B meson A = (e q ) B ~ 0.54 GeV [8]. 

Now we use the Gaussian approximation for the wave functions (13). Then shifting 
the integration variable p in (14) by — E £ ^ M A max , we can factor out the A dependence 
of the meson wave function overlap in form factor / + . The result can be written in the 
form 

/ + (0) = -F + (ALJexp(-CALJ, (17) 
where |A maa; | is given by (15) and 

CA 2 2A 2 AL, _A 2 ( M B -M n \ 2 

^ max (Pl + P^ + M^ f3 B '\M B + Mj ' 1 ] 

here 77 = „ 2 and A is equal to the mean value of light quark energy between B and ir 
meson states: 

A = (e q ) = —-^—e z K 1 (z) « 0.53 GeV, (19) 



2 

where z — — - 



Substituting the Gaussian wave functions (13) in the expression for the form factor 
(14) and taking into account (16)— (18), we get up to the first order in l/m b expansion: 



(3% J V m b \ f3 2 B n 3\e q + m q / \ + 

(-. ~~ rTj / (A + M w + 3m q ) + A77/2 (1 + — ) 

3(v/5 + 2) \{e q + m q ) 2 / V V \ \ ^ + 2/ 

' \(A + M w + 3m ff )-3 + i/ T7 T V) (20) 



,e q + m q / '' 3\e q (e q + m q ) 

\ / \ 3/2 

where N = ( pi^p J — \^V) * s due to the normalization of Gaussian wave functions 



in (21); e q = y p 2 + m 2 + A 2 r/ 2 , i. e. the energy of light quarks in final meson acquires 
additional contribution from the recoil momentum. 
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4 Results and discussion 



Using the parameters of Gaussian wave functions (13) in the expression (20) for the B — > 7r 
transition form factor /+(0) and in the similar expressions [15] for the B — > p transition 
form factors A^O), A 2 (0) and V(0) we get 



/f^(0) = 0.21 ±0.02 V B_ v(0) = 0.29 ± 0.03 
A^ p [o) = 0.26 ± 0.03 A^"(0) = 0.30 ± 0.03. 



(21) 



The theoretical uncertainty in (21) result mostly from the approximation of the wave 
functions by Gaussians (13) and does not exceed 10% of form factor values. We have also 
calculated the second order terms in l/m b expansion of form factors [15] and found them 
to be small (less than 5% of form factor values). 

We compare our results (21) for the form factors of B — > 7r(p)ez/ decays with the 
predictions of quark models [16,17] and QCD sym rules [18,19] in Table 1. There is 
an agreement between our value of f+~^ n (0) and QCD sum rule predictions, while our 
B — > peu form factors are approximately 1.5 times less than QCD sum rule results. 

To calculate the B — > n(p) semileptonic decay rates it is necessary to determine the 
g 2 -dependence of the form factors. Analysing the expressions (14), (17), (18), (20) for 
the form factor /+ and similar expressions [15] for Ai, A 2 and V, we find that the q 2 - 
dependence of these form factors near q 2 = is given by 

/+( ' 2) = ««< A -'' (21 > 



2JM B M I 

M^-^^-^ + ^HM^l (22) 



A2 ^ = In^tWMKn), (23) 

2JM B M p 



V ^ = (24) 



Mi+M z 



where w = 2M b mJ {p) i F+(*Lx) is defined b y ( 17 ) and A2(AL). V(A^J are 
defined similarly. We have introduced the function 

^Wn) exp h^nj' (25) 

which in the limit of infinitely heavy quarks in the initial and final mesons coinsides with 
the Isgur-Wise function of our model [8]. In this limit eqs. (21)-(24) reproduce the leading 
order prediction of HQET [9]. 

It is important to note that the form factor A\ in (22) has a different (^-dependence 
than the other form factors (21), (23), (24). In the quark models it is usually assumed 
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the pole [16] or exponential [17] g 2 -behaviour for all form factors. However, the recent 
QCD sum rule analysis indicate that the form factor A\ has g 2 -dependence different from 
other form factors [18,19]. In [19] it even decreases with the increasing q 2 as 

, , / q 2 \ , , , 2M B M a , \ t / \ / \ 

Mf) ^{i-^mo^ {MB ° M y i + «)M0). (26) 

Such behaviour corresponds to replacing £(u>) in (22) by £(w max ). 

We have calculated the decay rates of B — > ir(p)ev using our form factor values at 
q 2 = and the g 2 -dependence (21)-(25) in the whole kinematical region (model A). 
We have also used the usial pole dependence for form factors f + (q 2 ), A 2 (q 2 ), V(q 2 ) and 
Ai(q 2 ) = (m^+m P ) 2 (1 + w ) 1^2/^2 (model B), which corresponds to replacing the function 

£(w) (25) by the pole form factor. The results are presented in Table 2 in comparison 
with the quark model [16,17] and QCD sum rule [18,19] predictions. We see that our 
results for the above mentioned models A and B of form factor g 2 -dependence coincide 
within errors. The ratio of the rates T(B — > peu)/T(B — > neis) is considerably reduced 
in our model compared to the BSW [16] and ISGW [17] models, with the simple pole 
and exponential g 2 -behaviour of all form factors. Meanwhile our prediction for this ratio 
is in agreement with QCD sum rule results [18,19]. The absolute values of the rates 
F(B — > -neu) and T(B — > peu) in our model are close to those from QCD sum rules 
[19]. The predictions for the rates with longitudinally and transversely polarized p meson 
differ considerably in these approaches. This is mainly due to different g 2 -behaviour of A 1 
(see (22), (26) or usual pole dominance model [16]). Thus the measurement of the ratios 
F(B — > pev)jY{B — > neu) and Tl/Tt should provide the test of (^-dependence of A\ and 
may discriminate between these approaches. 

The differential decay spectra for B — > ir(p) semileptonic transitions, where x = 
u and Ei is lepton energy, are presented in Fig. 3 (see also [23]). 

We can use our results for V and A\ to test the HQET relation [20] between the form 
factors of the semileptonic and rare radiative decays of B mesons. Isgur and Wise [20] 
have shown that in the limit of infinitely heavy 6-quark mass an exact relation connects 
the form factors V and A 1 with the rare radiative decay B — > pq form factor Fi defined 
by: 

(p(p P , e)\uia^ u q u P R b\B(p B )) = ie^ Ta e* 1 'psPpF^q 2 ) 

+ [e;(M 2 B - M 2 ) - (e*q)(p B +p p ),}G 2 (q 2 ). (27) 
This relation is valid for q 2 values sufficiently close to q 2 nax = (M B — M p ) 2 and reads: 

2 l±Ml-Ml m + M B + M 
1W ' 2M B M B + M p 2M b 

It has been argued in [21,22,19], that in these processes the soft contributions dominate 
over the hard perturbative ones, and thus the Isgur- Wise relations (28) could be extended 
to the whole range of q 2 . In [10] we developed l/m& expansion for the rare radiative decay 
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form factor -Fi(O) using the same ideas as in the present discussion of semileptonic decays. 
It was shown that Isgur-Wise relation (28) is satisfied in our model at leading order of 
1/rrib expansion. The found value of the form factor of rare radiative decay B — > p7 up 
to the second order in l/m& expansion is [10] 

F?^ p (0) = 0.26 ± 0.03. (29) 

Using (28) and the values of form factors (21) we find 

F?^ p (0) = 0.27 ±0.03, (30) 

which is in accord with (29). Thus we conclude that l/m b corrections do not break the 
Isgur-Wise relation (28) in our model. 

5 Conclusions 

We have investigated the semileptonic decays of B mesons into light mesons. The recoil 
momentum of final 7r(p) meson is large compared to the ir(p) mass almost in the whole 
kinematical range. This requires the completely relativistic treatment of these decays. 
On the other hand, the presence of large recoil momentum, which for q 2 = is of order 
of m b /2, allows for the l/m b expansion of weak decay matrix element at this point. 
The contributions to this expantion come both from heavy 5-quark mass and large recoil 
momentum of light final meson. 

Using the quasipotential approach in quantum field theory and the relativistic quark 
model, we have performed the l/m& expansion of the semileptonic decay form factors at 
q 2 = up to the first order. We have determined the g 2 -dependence of the form factors 
near q 2 = 0. It has been found that the axial form factor A\ has a g 2 -behaviour different 
from other form factors (see (21)-(24)). This is in agreement with recent QCD sum rule 
results [18,19]. The ratios Y{B — > peu)/T(B — > ixev) and T L /T T are very sensetive to the 
(^-dependence of A 1 , and thus their experimental measurement may discriminate between 
different approaches. 

We have considered the relation between semileptonic decay form factors and rare 
radiative decay form factor [20], obtained in the limit of infinitely heavy 6-quark. It has 
been found that in our model l/m b corrections do not break this relation. 
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Table 1 Semileptonic B — > tx and B — > p decay form factors. 



Ref. 


/£-(0) 


A?^ p (0) 


A^ p (0) 


V B ^(0) 


our results 


0.21 ±0.02 


0.26 ±0.03 


0.30 ±0.03 


0.29 ±0.03 


16 


0.33 


0.28 


0.28 


0.33 


17 


0.09 


0.05 


0.02 


0.27 


18 


0.26 ±0.02 


0.5 ±0.1 


0.4 ±0.2 


0.6 ±0.2 


19 


0.23 ±0.02 


0.38 ±0.04 


0.45 ±0.05 


0.45 ±0.05 



Table 2 Semileptonic decay rates T(B — > neu), T(B — > pez/) (x|Kb| 2 x 10 12 s x ) and the 
ratio of the rates for longitudinally and transversely polarized p meson. 



Ref. 


F(B -> irev) 


F(B -> pei/) 


r L /r T 


our results 








model A 


3.1 ±0.6 


5.7 ± 1.2 


0.6 ±0.3 


model B 


3.0 ±0.6 


5.2 ± 1.2 


0.6 ±0.3 


16 


7.4 


26 


1.34 


17 


2.1 


8.3 


0.75 


18 


5.1 ± 1.1 


12 ±4 


0.06 ±0.02 


19 


3.6 ±0.6 


5.1 ± 1.0 


0.13 ±0.08 



Figure captions 

Fig. 1 Lowest order vertex function 

Fig. 2 Vertex function with the account of the quark interaction. Dashed line corresponds 
to the effective potential (4). Bold line denotes the negative-energy part of the quark 
propagator. 

Fig. 3 The differential decay spectra (model A) for B — > n(p) semileptonic transi- 
tions, where x = and E\ is lepton energy. Absolute rates ^ can be obtained using 
T(B — > Tr(p)eu) from Table 2. 
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This figure "figl-l.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9505321v2 



This figure "figl-2.png" is available in "png" format from: 



http://arXiv.org/ps/hep-ph/9505321v2 
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